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Abstract 

We discuss the sharp interface limit of a diffuse interface model for a two- 
phase flow of two partly miscible viscous Newtonian fluids of different densities, 
when a certain parameter e > related to the interface thickness tends to 
zero. In the case that the mobility stays positive or tends to zero slower than 
linearly in e we will prove that weak solutions tend to varifold solutions of 
a corresponding sharp interface model. But, if the mobility tends to zero 
faster than e 3 we will show that certain radially symmetric solutions tend to 
functions, which will not satisfy the Young-Laplace law at the interface in the 
limit. 
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1 Introduction 

The present contribution is devoted to the study of the relations between so-called 
diffuse and sharp interface models for the flow of two viscous incompressible Newto- 
nian fluids. Such two-phase flows play a fundamental role in many fluid dynamical 
applications in physics, chemistry, biology, and the engineering sciences. There are 
two basic types of models namely the (classical) sharp interface models, where the 
interface T(t) between the fluids is modeled as a (sufficiently smooth) surface and 
so-called diffuse interface models, where the "sharp" interface T(t) is replaced by 
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1 INTRODUCTION 



an interfacial region, where a suitable order parameter (e.g. the difference of volume 
fractions) varies smoothly, but with a large gradient between two distinguished values 
(e.g. ±1 for the difference of volume fractions). Then the natural question arises how 
diffuse and sharp interface models are related if a suitable parameter e > 0, which is 
related to the width of the diffuse interface, tends to zero. There are several results 
on this question, which are based on formally matched asymptotics calculations. But 
so far there are very few mathematically rigorous convergence results. 

More precisely, we study throughout the paper the sharp interface limit of the 
following diffuse interface model: 



f>0,v+ ( pv + 



• Vv - div(z/(c)£>v) + Vp = -e div(o(c) Vc <g> Vc) in Q, (1.1) 

clivv = inQ, (1.2) 

d t c + v ■ Vc = div(m e (c)Vp) inQ, (1.3) 

li = e^f'icj-eAc in Q, (1.4) 

v|on = on S, (1.5) 

nac • Vc| 9 q = n dU ■ Vp\ 9 n = on S, (1.6) 

(v,c)|t= = (v ,co) in Q, (1.7) 



where Q = Q x (0, oo),S = dQ x (0,oo), Q C ]R n is a suitable domain, and J = 
— m e (c)V/i. Here c = C2 — c\ is the volume fraction difference of the fluids, p = p(c) 
is the density of the fluid mixture, depending explicitly on c through p(c) = P2 ~ pi c — 
pl + p2 ; where pj is the specific density of fluid j = 1,2, and / is a suitable "double- 
well potential" e.g. /(c) = |(1 — c 2 ) 2 . Precise assumptions will be made below. 
Moreover, e > is a small parameter related to the interface thickness, p is the 
so-called chemical potential, m £ (c) > a mobility coefficient related to the strength 
of diffusion in the mixture and a(c) is a coefficient in front of the |Vc| 2 -term in the 
free energy of the system. Finally, iiqq denotes the exterior normal of dQ. The model 
was derived by A., Garcke, and Griin [5]. In the case p(c) = const, it coincides with 
the so-called "Model H" in Hohenberg and Halperin [14], cf. also Gurtin et al. [13] . 
Existence of weak solutions for this system in the case of a bounded, sufficiently 
smooth domain Q and for a suitable class of singular free energy densities / was 
proved by A., Depner, and Garcke [I]. We refer to the latter article for further 
references concerning analytic results for this diffuse interface model in the case 
p(c) = const, and related models. 

In [5] the sharp interface limit e — > was discussed with the method of formally 
matched asymptotics. It was shown that for the scaling m £ (c) = rhe a with a = 0, 1, 
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m > 0, solutions of the system (1.1)-( L5| converges to solutions of 

p ± d t v + (p ± v + ?^J)- Vv-divT^v,?) = mQ ± (t),t>0, (1.8) 

divv = mQ ± {t),t>0, (1.9) 

m A/2 = infi ± (t),t>0, (1.10) 

-n- [T(v,p)] = aHn onr(t),t>0, (1.11) 

V - n ■ v|r(t) = -[^n ■ V/i] on r(t),t > 0, (1.12) 

^| r(t) = o-ff onr(t),t>0, (1.13) 

with J = — m V/i. Here n denotes the unit normal of T(t) that points inside Q + (t) 
and V and H the normal velocity and scalar mean curvature of T(t) with respect to 
n. Moreover, by [•] we denote the jump of a quantity across the interface in direction 
of n, i.e., [f](x) = \im h ^ (f(x + hn) — f(x — hn)) for x G T(t). Furthermore, a is 
a surface tension coefficient determined uniquely by / and m = fh if a = and 
m = if a = 1 is a mobility constant. Implicitly it is assumed that v,/x do not 
jump across T(t), i.e., 

[v] = [/i] = onr(t),t>0. 
In the following we close the system with the boundary and initial conditions 

v\ d n = ondQ,t>0, (1.14) 

nan ■ m e (c)V/i|an = on 90, t > 0, (1-15) 

O+(0)=O+, (1.16) 

v| t=0 = v in Q, (1.17) 

where v ,Oq are given initial data satisfying dfl^ fl dQ = 0. Equations (1.8 )-( 1.9 
describe the conservation of linear momentum and mass in both fluids, and (1.11 



is the balance of forces at the boundary. The equations for v are complemented by 
the non-slip condition (1.14) at the boundary of Q. The conditions (1.10), (1.15) 



describe together with (1.12) a continuity equation for the mass of the phases, and 
(1.13) relates the chemical potential \i to the L 2 -gradient of the surface area, which 



is given by the mean curvature of the interface. 

We note that in the case a = 1, i.e., m = 0, (1.12) describes the usual kine- 



matic condition that the interface is transported by the flow of the surrounding 



fluids and ( 1.8 )-( 1.17) reduces to the classical model of a two-phase Navier-Stokes 



flow. Existence of strong solutions locally in time was first proved by Denisova and 
Solonnikov jTTJ. We refer to Priiss and Simonett [19] and Kohne et al. [H] for more 
recent results and further references. Existence of generalized solutions globally in 
times was shown by Plotnikov [18] and A. [1, 2]. On the other hand, if a = 0, mo > 0, 



respectively, the equations (1.10), (1.13), (1.15) are a variant of the Mullins-Sekerka 



flow of a family of interfaces with an additional convection term n-v|r(t)- In the case 
Pi — Pi existence of weak solutions for large times and general initial data was proved 
by A. and Roger [6] and existence of strong solutions locally in time and stability of 
spherical droplets was proved by A. and Wilke [8]. 
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2 NOTATION AND PRELIMINARIES 



In the following we address the following question: Under which assumptions 
on the behavior of m £ (c) as e — » do weak solutions of ( 1.1 )-( 1.7) converge to 
weak/generalized solutions of (1.8)- (1.17)? In this paper we provide a partial answer 
to that question. If one assumes e.g. m £ (c) = rhe a , the results in the following will 
show that convergence holds true in the case a G [0, 1). More precisely, we will show 
that weak solutions of ( 1.1 )-( 1.7) converge to so-called varifold solutions of (1.8)- 
(1.17), which are defined in the spirit of Chen [TU]. But in the case a G (3, oo) we 
will construct radially symmetric solutions of (1.1 )-( 1.4) in the domain Q = {x G ffi. : 
1 < \x\ < M} with suitable inflow and outflow boundary conditions, which do not 
converge to a solution of ( 1.8 )-( 1.13). In particular, the pressure p in the limit e — > 
satisfies 

[p] = an(t)H onT(t) = dB m (0), 

where R(t), n(t) — >t->oc 00 an d v is independent of t and smooth in Q. This shows 
that the Young-Laplace law (1.11) is not satisfied. We note that these results are 
consistent with the numerical studies of Jacqmin, where a scaling of the mobility as 
m e (c) = rhe a with a G [1,2) was proposed and considered. 

The structure of the article is as follows: First we introduce some notation and 
preliminary results in Section [2j Then we prove our main result on convergence of 
weak solutions of ( 1.1 )-( 1.7) to varifold solutions of (1.8 )-( 1.17) in the case that the 
mobility m B (c) tends to zero as e — > slower than linearly in Sectio n [3] Finally, in 
Section |4j we consider certain radially symmetric solutions of ( 1.1 )-( 1.7) and show 
that these do not converge to a solution of ( 1.8 )-( 1.13) if the mobility tends to zero 
too fast as s — » 0. 



2 Notation and Preliminaries 

Let U C R d be open. Then M(U;R N ), N > 1, denotes the space of all finite Re- 
valued Radon measures on U. By the Riesz representation theorem M.(U;M. N ) = 
C (U;R N )', cf. e.g. Theorem 1.54]. Moreover, M{U) := M(U,R). Given 
A G M.(U;R ) we denote by |A| the total variation measure defined by 

{oo oo 
^2 \K A k)\ ■ A k G B{U) pairwise disjoint, A = (J A k 
k=0 k=0 

for every A G B(U), where B(U) denotes the cr-algebra of Borel sets of U. Moreover, 
U — > R N denotes the Radon-Nikodym derivative of A with respect to |A|. The 
restriction of a measure fj, to a /i-measurable set A is denoted by (/j i \_A)(B) = n(Ar\B). 
Furthermore, the s-dimensional Hausdorff measure on R d , < s < d, is denoted by 
T-L s . Recall that 

BV(U) = {/ G L\U) : V/ G M(U; R d )} 
\\f\\BV(U) = Wfh^U) + ||V/||M(f7;R<i), 
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where V/ denotes the distributional derivative. Moreover, BV(U;{0, 1}) denotes 
the set of all X G BV(U) such that X(x) G {0, 1} for almost all x G U. 

A set E C [/ is said to have finite perimeter in Z7 if Afg G BV(U). By the 
structure theorem of sets of finite perimeter |VA^| = 'H d ~ 1 [d*E, where d*E is the 
so-called reduced boundary of E and for all (p G C (U, M. d ) 

-(VX E ,<p) = / divpdx = - I y-n E dU d ~ l , 

JE Jd*E 

where n E (x) = iy^i , cf. e.g. [9]. Note that, if E is a domain with (^-boundary, 
then d*E = dE and coincides with the interior unit normal. 

As usual the space of smooth and compactly supported functions in an open set 
U is denoted by Cq°(U). Moreover, C°°(U) denotes the set of all smooth functions 
f'.U—t C such that all derivatives have continuous extensions on U . For < T < oo, 
we denote by Lf oc ([0,T);X), 1 < p < oo, the space of all strongly measurable 
/: (0,T) -> X such that / G L p (0,T';X) for all < T < T. Here L P (M) and 
L P (M; X) denote the standard Lebesgue spaces for scalar and X-valued functions, 
respectively. Furthermore, C^° CT (f2) = {cp G C£°(£l) d : div cp = 0} and 

If F = X' is a dual space and Q C R N is open, then L^(Q; F) denotes the space of 
all functions z/: Q — > Y that are weakly-* measurable and essentially bounded, i.e., 

x ^ (u x ,F(x, -))x',x 
is measurable for each F G L l (Q\ X) and 

M\l~(Q;Y) ■= eSS SUp x£ Q||l/ x ||y < OO. 

Moreover, we note that there is a separable Banach space X such that X' = BV(Q), 
cf. 0. As a consequence [12] we obtain that L%{0,T;BV(Q)) = (L\0,T;X))* and 
that uniformly bounded sets in L^,(0,T; BV(Q)) are weakly *-precompact. 



3 Sharp Interface Limit 



In this section we discuss the relation between ( 1.8 )-( 1.17) and its diffuse interface 
analogue ( 1.1 )-( 1.7). 

Assumption 3.1 We assume that the domain Q C M d , d = 2,3 is bounded and 
smooth. Furthermore, we assume that there exist constants cq,Co > such that 



f G C 3 (R), /(c) > 0, f(c) =0 if and only if c 
|c| > 1 — Co for some constant p > 3 



-1,1, and f"(c) > c \c\ p - 2 if 
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3 SHARP INTERFACE LIMIT 



p,a,u G C 1 (IR) with c < p,a,v < C and 

P2 + Pi , Pi - Pi 



P(c) 



for c G [-1, 1] 



• m e ,m G C 1 (IR) ; < m e ,m < Co, m £ — > £ ->o ^ n C 1 ^); a nd either mo > Co 
or m = 0. J/mo = ; then m £ > fn £ for constants m £ > with m £ — t- £ ^ 0. 

The stronger assumption p > 3 (compared to p > 2 in [10]) is needed here 
for the uniform estimate of v e • Vc £ = div(v e c e ) in L 2 (0, T; A possible 

choice for the homogeneous free energy density is f(s) = (s 2 — l) 2 . Moreover, let 
a = y/ f(s)/2 ds and A(s) = f Q s a/ a(r) dr. 

Now, let us consider the energy identities corresponding to our two systems. We 
recall that every sufficiently smooth solution of the Navier-Stokes/Mullins-Sekerka 
system ( 1.8 )-( 1. 17) satisfies 



d 1 
dt2 



[ p{c) \v\ 2 dx + a < ^-H d - 1 {Y) = - [ v{c)\Dv\ 2 dx- [ m (c)| V/i| 2 dx, (3.1) 
Jn dt Jq J n 

where c(t, x) = —l+2xn+<t)( x )- On the other hand, every sufficiently smooth solution 



of ( 1.1 )-( 1.7) satisfies 



d 1 
~dt2 



J p{c) \\\ 2 dx + ^£ £ 



- / v(c)\Dv\ 2 dx 



m £ (c)\Vp\ 2 dx, (3.2) 



where 



SJc) 



JVA(c)\ 2 | f(c) 



dx 



2 : 

is the free energy. Moreover, by Modica and Mortola [T7j or Modica [16] . for A(c) = c, 
we have 

£ £ — 'P w.r.t. L 1 -r-convergence, 

where 



V{u) 



al-i d 1 (d*E) if u = —1 + 2\e and E has finite perimeter, 
-oo else. 



Here, d*E denotes the reduced boundary. Note that d*E = BE if E is a sufficiently 
regular domain. Therefore, we see that the energy identity (3.1) is formally identical 
to the sharp interface limit of the energy identity (3.2) of the diffuse interface model 

We will now adapt the arguments of Chen [TU], see also A. and Roger |6J, to show 
that, as e — > 0, solutions of the diffuse interface model (1.1 )-( 1 .7 ) converge to varifold 
solutions of the system ( 1.8 )-( 1.17). Let Q = Q x (0, oo) and Gd-i '■= S d ~ l / ~ where 
uo ~ v\ for uq, V\ G S d ~ x iff v§ = ±v\ and S d ~ l is the unit sphere in IR d . 
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Definition 3.2 Let v G L^(Q) and E C Q be a set of finite perimeter. Then 



(y,E,fi,V) if m > and (v,E,V) else is called a varifold solution of (l.S)-(l.ll) 
with initial values (vo, Eq) if the following conditions are satisfied: 

1. v G L 2 ((0,oo);i/ 1 (fi) d )nL oo ((0,oo);L2(fi)); p e Lf oc ([0, oo); , V/i G 
L 2 ((0,oo);L 2 (fi) d ) ifm > 0. 

2. E = \J t>0 E t x {t} is a measurable subset of Q x [0, oo) swcn i/iai G 
CQC^oo/fL 1 ^)) nL~.((0,oo);BV(n)) and |£ t | = |£7„| /or a// 1 > 0. 

3. V is a Radon measure on Q x G^-i x (0, oo) such that V = V t dt where V 1 is a 
Radon measure on Q x G^-i / or almost all t G (0, oo), i.e., a general varifold 
in Q. Moreover, for almost all t G (0, oo) V 1 has the representation 

f ^(x 1 p)dV t (x,p) = V fb\{x)^{x,p\{x))d\\x) (3.3) 
JnxGd-! i=1 Jn 

for all ip G C(Q x Cr^-i)- i/ere, /or almost all t G (0, oo), A* is a Radon 
measure on Q, and the A* -measurable functions b\,p\ are R- and G^-i-valued, 
respectively, such that 

d d 
< 6| < 1, > 1, Y,P\®P\ = I A< "«- e - 

i=l i=i 

and 

A* - 2a' 

^. For c := — 1 + 2x^, J := — m (c)V/z i/ m > and J := e/se as well as 
J := §^(c)J we nave 

y ( ~ P( c ) v • - v ® (p(c)v + J) : Vy> + z/(c)L>v : D<p) d(x, i) 

- / p(c| t=0 ) v ■ ip\ t=Q dx = - (5V\ if) dt (3.4) 
Jn Jo 

for all if G Cq°([0, oo); C^(fi)) and 

2 / + div(^v) d(x, t) + / J -Vil>d(x,t)+ / ^|t=odx = (3.5) 

J E Jq Je 

for all ip G C °°([0, oo) x H). Here 

(8V\<p):= f_ (I-p®p) : V^d(x,p) /or a// G C°°(n ; M d ). 
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3 SHARP INTERFACE LIMIT 



Furthermore, if m > we have 

2 / div(//77) dx = (5V\ 77) (3.6) 

JEt 



for all r) £ Cq(Q; IR d ) and almost all t £ (0, oo) 
5. Finally, for almost all < s < t < oo 



p{c{t))\v(t)\ 2 dx + \ t (n)+ I I {v(c)\Dv\ 2 - J V/i) d(x,r) 

< s / p(c(s))|v( S )| 2 da; + A s (n). (3.7) 



/ / (z/(c)|Dv| 2 - J • V/i) d(av 
is in 



n 



2 

We define the free energy density by 

|V^(c)| 2 f(c) 

2 £ 

In [3] the existence of global weak solutions is shown for a class of singular free 

energies. We note that this proof can be easily carried over to the present situation 
with only minor modifications and even some simplifications since / is non-singular. 
Throughout this paper we will use the definition of weak solutions in [3J. By this 
definition we have 

v £ £ BC U ([Q, oo); Ll (O)) n L 2 ((0, oo); H\Q) d ), 

c £ £ BC U ([Q, oo); H\n)) n L 2 OC ([0, oo); tf 2 (fi)), /(c e ) £ L 2 OC ([0, oo); L 2 (fi)), 
// £ L 2 OC ([0, oo); L 2 (fi)), V/i £ L 2 ([0, oo); L W), 

and 



/ -p(c £ )v £ ■ <9 t c^ - v e ® (p(c £ )v £ + J B ) : Vc^> + u(c £ )D\ £ : Dp d(x, t) 
JQ 

- \ p(c , £ )v 0!£ ■ cp\ t=0 dx = / ea(c £ )Vc £ <g> Vc £ : Vipd(x,t) (3.8) 
in iQ 

for J £ := -m £ (c £ )V/i £ , J £ := f (c e )J e , and all <p £ C o °°([0, oo); C£(fl)), as well as 

c £ (<9tV> + div(^v £ )) d(x,t) + / c 0i£ ^| 4=0 (ia; = / m E (c e )V/i E -Vip d(x,t) (3.9) 

in ig 

for all i) £ C$°([0, oo) x U), and 

/i £ = — + ea'(c £ ) — ^ e div(a(c £ )Vc £ ) a.e. in Q, (3.10) 

nan • Vc £ = a.e. on (0, oo) x <9f2. (3.11) 
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Moreover, we have 



P(Ce(t)) J V£(t)|2 cix + g £ (c £ (t))+ f [ v(c £ )\Dv £ \ 2 dxdr 
u * Js Jn 

+ [ [ m £ (c £ )|V/i £ | 2 d(x,t) < / PM^Stl^ll dx + S e (c e (s)) (3.12) 
is Jn Jn A 

for all t > s and almost every s > including s = 0. 

Theorem 3.3 For a// e G (0, 1], Zei initial data (vo )E , co,e) G -^(fi) x -fZ^fi) fre gzi>en 
suc/i t/iat p| J n Co j£ dx = c G (—1, 1) and 

/* IZ|i^ rfx + £ £ ( COj£ ) < R (3.13) 



/or some R > 0. Furthermore, let (v e , c e ,// £ ) &e weaA; solutions of (l.l)-(l.l) in the 
interval [0, oo). Tnen tnere exzsis a sequence (ek)keN> converging to as k oo, 
such that the following assertions are true. 

1. There are v G L 2 ((0, oo); # 1 (fi) d ) n L°°((0, oo); L 2 (fi) d ), v G L 2 (fi) sncn that, 
as k — > oo ; 

v £fc -v m J L 2 ((0,oo);tf 1 (fi) d ), (3.14) 
v £fe ->v mL 2 oc ([0,oo);L 2 (fi)), (3.15) 
v 0i£fc -v mL 2 (fi). (3.16) 

J/m > 0, there exists a fi G £ 2 OC ([0, oo); /^(fi)) with V/x G L 2 ((0, oo); L 2 (fi) d ) 
and stzc/i inai 

/i £fe -// mLLaO.oo);^)). (3.17) 
There are measurable sets E C x [0, oo) and i?o C fi st/cn inai, as k — )■ oo, 

c £fc ->■ -1 + 2xe a.e. m fi x (0, oo) and m C£ c ([0, oo); L 2 (fi)) (3.18) 
c o,e fc — > — 1 + 2Xe a - e - i n fi- (3.19) 

In particular, we have XE\t=o — Xe in L 2 (Q). 

3. There exist Radon measures A and A^-, i, j = 1, . . . , d on fi x [0, oo) such that 
for every T > 0, i, j = 1, . . . , d, as k — ^ oo, 

e £k {c £k ) dx dt ^* X in M{U x [0,T]), (3.20) 
e k a(c £k )d Xi c £k d Xj c £k dxdt ^* Ay in -M(fi x [0,T]). (3.21) 
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3 SHARP INTERFACE LIMIT 



4- There exists a Radon measure V = V t dt on Q x G^-i x (0, oo) such that 
(v, E, fi, V) ifm > and (v, E, V) else is a varifold solution of (1.8)-(1.17) in 
the sense of Definition 3.2 with initial values (v , E ) and a = J\ \J f(s)/2ds. 
Furthermore, 

[ (SV t ,rj)dt= [ [ Vrj: (dXI-idXij)^) dt (3.22) 
Jo Jo Jn 

for all riE C^tt x [0,T];M d ). 

5. //v 0)£fc —> v in Ll(Q) and £ £ (co, e ) ~^ 2cr|Vx_B |(^) as k ^ oo, then (3.7) holds 
for almost all t G (0, oo), s = 0, and A°(f2) replaced by 2cr|Vx£ l(^)- 



By (3.12) and the assumptions on the initial data we obtain 

p(c £ (t))\Mt)\ 2 



+ 



dx + S e (c e (t)) 



v(c £ )\Dv £ \ 2 + m £ (c e )|V fi e \ 2 dx dt < R 



(3.23) 



o ./n 



for all t > 0. 

From this estimate, Korn's inequality, and (3.13) we deduce that there exists a 
sequence \ as k — > oo such that (3. 14), (3. 16), (3. 17), (3. 20), and (3.21) hold. 
Using the assumptions on /, we further deduce that 

/ \c £ {t)\ p dx < C{l + R), (3.24) 
Jn 

(\c e (t)\ - Ifdx < CeR (3.25) 



for all t > 0. In particular, for (3.25) we used that f(c) > C(\c\ — l) 2 for all c G K and 
some constant C > which follows from the positivity of /"(±1) and the p-growth 
of / for large c. With the definitions (cf. [10J) 



W(c) 



2/(s) ds, where f(s) = min(/(s), 1 + \s\ 



and 



w £ (x,t) = W(c £ (x,t)), 
the functions w £ are uniformly bounded in L°°((0, oo); BV(Q)) since 



\Vw e (x,t)\dx= / A/2/(c £ (x,t))|Vc e (x,t)|rfx < / e £ (c £ (x, t)) dx < R. (3.26) 
Jn Jn 

Moreover, note that by the assumptions on /, there exist constants C , Ci > such 
that for all Cq, c\ Gl 

Co|c - Cl | 2 < |^(co) - W( Cl )| < Ci|co - Cl |(l + [co| + | Cl |). (3.27) 

Here, for the first inequality we used again that f(s) > C(|s| — l) 2 for all s G R. 
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Lemma 3.4 There exists a constant C > such that 



\w 



+ \\Ce 



£|l C5([0,oo);L 1 (Q)) 11 eil CS([0,oo);L 2 (Q)) 



< c. 



Proof: The proof is a modification of PJjJ Proof of Lemma 3.2]. Therefore, we only 
give a brief presentation. For sufficiently small rj > 0, x G f2, and t > let 



uj(y)c e (x -r]y,t)dy, 



where u is a standard mollifying kernel and c £ is extended to a small neighborhood 
of Q as in [TUl Proof of Lemma 3.2]. Then, there exist constants C, C > such that 

II Vc£(*) Hi^n) < CV'IMt)!^ < CV 1 (3.28) 
\\4(t) - c £ (t)f L2m < CT,\\Vw e (t)\\»W < C'v (3.29) 



for all sufficiently small i] > 0, cf. [TOT Proof of Lemma 3.2]. From (3.28) and (3.9) 
we deduce that for all < r < t < oo such that \t — r| < 1 

(c e (a;,t) - c e (x, T))(c v £ (x,t) - c v £ (x,T))dx = 

- I / (m £ (c £ (x, s))Vfi E (x, s) - v £ (x, s)c £ (x, s)j ■ (Vc*(x, t) - Vc v £ (x, r) j d(x, s) 
<C(R)(t-r)* sup HVc^^ll^^C^-^t-r)!. (3.30) 

se(r,t) 

Here, we used the fact that for all r, t as above we have 



|m £ (c £ )V/i e - v e c e 



IL 2 (Qx(-r,t)) 



< C(i2), 



since the sequences (v e ) C £ ((0, oo); L 6 (Q)) and (c e ) C L°°((0, oo); L 3 (fi)) are 



bounded due to the assumptions d < 3 and p > 3. Now, combining (3.30), (3.29) 



and using Holder's and Young's inequality we conclude that for rj, r, and t as above 
we have 

\\Ce(t) - C e (r)\\ 2 LHn) <C(t7 + tT 1 |*-T|*). 

Choosing 77 = (£ — r) 2 for sufficiently small t — t we conclude the claim concerning 



c £ . Using (3.27), one derives the claim concerning w £ as in [TO] . 



Remark 3.5 It is possible to understand the proof of Lemma [3.4| from a more gen- 
eral point of view 

In particular 



time- derivative of (c 



From (3.9) and (3.23) we easily deduce that the distributional 
is uniformly bounded in L 2 ((0,oo) 



H 



(c e ) is uniformly bounded in C 1//2 ([0, oo), H x (fi)). On the other hand, the computa- 



1 /3 

tions leading to (3.29) show that (c e ) is uniformly bounded in L°°((0, oo); (ft)). 
This follows from B~^(Q) = (L 2 (fi), if 1 (r2)) 1 / 3i00 and the definition of the real inter- 
polation spaces with the aid of the .fT-method. By interpolation, we obtain uniform 
boundedness in C x ^([0, oo), L 2 (fi)). 
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The proof of the following lemma is literally the same as the proof of [TUl Lemma 3.3] . 

Lemma 3.6 There exists a subsequence (again denoted by Sk) and a measurable set 
E C Vt x [0, oo) such that, as k — > oo, 

i 

w £k — > 2o~Xe o-e. in x (0, oo) and in C^ c ([0, oo); 

i 

c Ek -»■ -1 + 2xe a.e. in Q x (0, oo) and in C£ c ([0, oo); L 2 (Q)) 
Moreover, X e G L~((0, oo); BV(Q)) H C*([0, oo); L 1 ^)) and /or all t > we have 

\E t \ = \E \ = ^\n\. 

Lemma 3.7 There exist constants C, Eo > such that 

\\Ve(t)\\m(n) < C (8 £ (c £ (t)) + ||V^ £ (*)|U 2(n) ) (3.31) 



for almost all t > and < e < sq. Using m £ > m £ we deduce from (3.31) and 



m £ 



I \\He{t) lliarn) dt ^ T ) /° r all < T < oo. (3.32) 
•/ii 

Proof: Let us suppress the time variable. Due to Poincare's inequality it suffices to 



control the average of /i £ . Equation (3.10) can be written in the form 



f\c e ) 



He = - eVa(c e ) AA(c E ). (3.33) 

Multiplying by t] ■ Vc £ for rj G C 1 (fi;IR d ), integrating over Q, and integrating by 
parts yields 

/ rj • Vc £ fi e dx = - / Vr7 : (e £ (c £ ) I — e V A{c £ ) <g> Vi(c E )) da; 

+ / e e {c £ )-n-n d adU d - 1 . (3.34) 
Jan 

Now we can proceed exactly as in the proof of JTOj Lemma 3.4]. ■ 

Lemma 3.8 There exists a subsequence (again denoted by Sk) such that, as k — )■ oo ; 

v £fc ^v znLL([0,oo);L^(fi)) 
v £jb (t) — y v(t) in L 2 (fi) for almost every t > 0. 

Furthermore, there exists a measurable, non-increasing function £{t), t > 0, such 
that for almost all t > 

£ £k (c £k (t)) £(t) and |Vx^|(fi) < ^£(t) < ^J2. (3.35) 
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Proof: Let us fix some T > and let P a : L 2 {Q) d — > L%(Q) denote the Helmholtz 
projection. In order to prove the claim concerning v £fc it suffices to show that for a 
subsequence we have 



PMOv*) Pa(p(c)v) in L'((o,r); (K(n) n # W)') (3.36) 

since then 

/ / P(ce k )\v ei .\ 3 ddt= [ [ P CT (p(c £ Jv £ J -v £k dxdt 
Jo Jn Jo Jn 

-^-kx> / / P CT (p(c)v) • v dx dt = / / p(c)|v| 2 dxdt, 
Jo Jn Jo Jn 

and from this convergence, the strong convergence of (c £fe ), and the strict positivity 
of p we easily deduce the claim, cf. [3j. But (3.36) follows from the Aubin-Lions 
lemma by noting that, firstly, 

L 2 a (9) (Ll(n) n H 1 (n) d y ^ (Ll(n) n w^°°(n)y 

and that, secondly, the distributional time-derivative of (P a (p(c £k )v £k )) is uniformly 
bounded in L 8 / 7 ((0,T); n W 1 ' 00 ^))')- This last bound follows by estimating 

each term in (Q. We have (appreviating L p ((0, T); by L p L q ) 

||p(c e )v e <g> v £ || L 2 L 3/2 < ||p(c e )v £ || i0 o L 2||v £ || L 2 L 6, 

||V £ ® J £ || L 8/7 L 4/3 < ||V £ ® J £ ||^ 4 L3/2 ||V £ g) J e ||^il 

< C||v £ ||5 / / 2 4 L6 ||v £ ||J / / i L 2||m(c £ )|V / u £ | 2 || L i L i, 

||Z/(C £ )L>V £ || L 2 L 2 < C||L)V £ || L 2 L 2, 

||ea(c £ )Vc £ g> Vc £ || L oo L i < C\\e\V A(c £ )\ 2 \\ L o. L i. 
Concerning the remaining claims we note that the total energies 

sr{t)-.= \\\vM\h m + UcM\ t>o, 

form a sequence of bounded, non-increasing functions and that v £fc (t) — >k-+oa v (^) 
for almost all t > in L 2 (Q). Now, we can proceed exactly as in the proof of [TtH 
Lemma 3.3]. ■ 

Finally, we define the discrepancy function by 

C(ce) := £ -\VA(c £ )\ 2 -- £ f(c £ ). 
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Theorem 3.9 For all sufficiently small 77 > there exists a constant C{rj) such that 
for all sufficiently small e > ( the maximal e may depend on rj) we have 



[ (e(c £ )) + d(x,t)<r] [ [ e E (c £ )d(x,t)+EC(r)) [ [ |/i e | 2 d(x, t). 
Jn Jo Jvl Jo Jn 



Combining this estimate with the assumption e/m £ — > e ^o and (3.32) we deduce 
that 

Inn / / (f (c e )) + d(x, t) = for all < T < 00. 



Jn 



Proof: The proof is based on the elliptic equation (3.33) which can be written in 
the form 



Let c± := A(±l), B(c) := c c -^- + ^+±^, and f\c) := f(A^ 1 (B(c)))/(B') 2 for c G K. 
Then / fulfills Assumption 3.1, and for c £ := B~ 1 (A(c £ )) we have 



— e Ac F 



Since the function a(A(c £ )) l l 2 (B') 1 is uniformly bounded, [TQl Theorem 3.6] yields 



'0 J si 
where 



(f (c £ )) + d(x,t) < v I / e £ (c £ )d(x,t) + eC(rj) / / \fM £ \'d(x,t) (3.37) 
Jn Jo Jn 



f(5 e ) 
~e £ (ce) 

This proves the claim. 



||VC £ | 2 - -J{c £ )=nCe)/{B 

l\V~c £ \ 2 + 1 -f(5 £ )=e%c £ )/(B') 



f\2 



f\2 



Using the previous statements, we can now easily finish the proof of Theorem 3.3 



by the arguments of [TDJ Section 3.5]. To be more precise, item 1 follows from (3.13) 



and Lemmas 3.7 and |3.8[ Item 2 follows from Lemma 3.6 and the energy inequality 
(3.23). Item 3 follows from (3.23) as well. Furthermore, we note that A = X t dt for 



Radon measures A on Q since 



X(A x I) < A (ft x J) = lim / S £k (r) dr < \I\R 
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for any measurable A C Q, I C [0, oo). Similarly, we also get A*(fi) = S(t) for almost 
all t G (0, oo) due to (3.35). From (3.12) we deduce that 

A'(H)= lim £ £fc (c £fc (t)) 



Ai— >oo 



< 



-Hiriinf / / (z/(c £ J|£)v £fc | 2 + m £fc (c £ J |V^ £ J 2 ) d(x, r) 
+ lim ( ^e k (c £k (s)) + - / p(c £k (s))\v £k (s)\ 2 dx - - / 
- / / (v{c)\Dv\ 2 dx - J V/i) rf(x,r) + A S (H) 



P(c £fc W)|v £fc (t)| 2 rfx 



+ 



p(c(s))|v(s)| 2 <ix 



p(c(t))\v(t)\ 2 dx 



for almost all < s < t < oo where c := — 1 + 2%£. This is (3.7). Item 5 follows sim- 
ilarly. We can proceed as in [101 Section 3.5] to c onst ruct the varifold V. Therefore, 
we only give a sketch. We deduce from Theorem 3.9 that for all T] ,r] l G C 1 (fi;lR d ) 
and all < T < oo 

/ jvo^li ■ (dXij) < [ f_\Vo\\m\ dX. 
Jo Jn Jo Jn 

This proves the existence of A-measurable IR-valued, non-negative functions 7, and 
A-measurable unit vector fields % — 1, . . . , d, such that 

d d d 

(Xij) =2_,li l, i® an d /^7t < lj Uj® Uj = I A-a.e. 



1=1 



i=l 



We denote the equivalence class of Ui(x,t) in Gd-\ by p'(x), define the functions b\ 
by 



1 

b\{x) := 7i (x, t) + — (l - *)) 

i=i 



and define the varifold V as in (3.3). Then item 3 in Definition 3.2 follows taking 
into account (3.35). Furthermore, in the case mo > we infer from (3.34) that 

/ 2 XEt tiv(m)dx= [ Vr):(dXI-(dX tJ )* J=1 )= [ Vry : V b\(l - p\ ® p*) dX 
Jn Jn Jn i=1 

= (5V t ,r ] ) 

for all 77 G CQ(fi;lR d ) and almost all t G (0,oo). This is (3.6). Furthermore, these 
calculations prove (3.22). Similarly, (3.4) and (3.5) follow from (3.8) and (3.9), 
respectively, where one uses that 

/ ea(c £ ) Vc £ ® Vc £ : Vc^ d(x, t) = / ip ■ Vc £ /i £ d(x, t) ^ £ ^ (6V*, ip) 
Jo Jo 
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for all tp G C°°([0, oo); C£°(fi)). This proves item 4 in Definition [3~2| Finally, item 2 
in Definition |3.2 follows from Lemma 3.6 This concludes the proof of Theorem 3.3 



In the radially symmetric case we can prove a stronger statement concerning the 
discrepancy measure. 

Theorem 3.10 Let Q = Bx(Q), and assume that the solutions (v e , c £ , fi E ) are radially 
symmetric. Assume, furthermore, that A(c) = c for all c G 1R, and that the constants 



m £ in the Assumptions 3.1 satisfy 



£ d-i j me 



£-5-0 



0. 



(3.38) 



Then, for all T > 0, we have 



T 



lim / / \e(c £ )\d(x,t) = 0. 



For the proof we need the following result from [mi Lemma 4.4]. 

Lemma 3.11 There exist positive constants C and r/ such that for every r\ e [0, 770] 
e e (0, 1], and every (u £ ,v £ ) G H 2 (Q) x L 2 (Q) such that 



-sAu £ + e- l f(u^ 



nan • Vw £ 



\dn 







we have 



{xeC:M e |>l-»7} 



\u £ )+E-\f\u £ )Y) 



e\Vu e \ 2 dx + C e / \v e \ 2 dx 

{x£tl:\ut\<l-r)} JO. 



(3.39) 



Proof of Theorem 3.10: We can show exactly like in [10, Proof of Theorem 5.1] 
that there exists a constant C > such that for almost alH > we have 



(3.40) 
(3.41) 



e £ (c £ (t)) dx < C5M £ (t) for all 5 G (0, 1), 
|f (c e (r,t)) + fi E (r,t)c e (r,t)\ < Cr^M^t) for all r G (0,1). 
Here, we use the notation r = |sc| and 

M e (t) := 1 + \\fMe(t)\\m m +e\\fM £ (t)\\ 2 Hl(Q) . 



From (3.41) we deduce that for small 5, rj > 



\C(c £ (t))\dx< e £ (c £ (t))dx+ \fi e (t)\(l - r}) dx 

Q JB s U{\c e (t)\>l-ri} J nn{r>6,\c £ (t)\<l-ri} 



+ CM e (t) I 

/nn{r><5,|c £ (t)|<l-r ? } 



l-d 



dx. 
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Using (3.40) and (3.39), the first integral on the right hand side may be estimated 
by 

C5M%t) + C' V £ £ (t) + C'e\\fi £ (t)\\h m - 



By (3.25), the second integral is dominated by 

\ LHn) \{\c £ (t)\>l-v}\ 1/2 <C"( V )M^ty/ 2 . 



Finally, using (3.39) again, the third integral is smaller than C {rf)M e {t)5 1 d s. Sum- 
ming up, we have 



\e(ce(t))\dx < C'r]£ e (t) + C'e\\^(t)\\l* m + C"(r])M%t)(e^ 2 + S^e + 5). 



Integrating this estimate from to T and choosing rj small, the first term on the 
right hand side gets arbitrarily small. Choosing then S = e 1 ^ 2d ^ and e small the 
other two terms get arbitrarily small, too. While this is obvious for second term, 
concerning the third term we remark that it takes the form 



C"(rj) / M £ (t) dt {e l/2 + s 1/( - 2d - 2) ) = o(l) 
Jo 



as £ — > 



due to (3.38). 



4 Nonconvergence 



In this section we show that solutions of ( 1. 1 )-( 1.4 ) do not converge in general to solu- 
tions of ( 1.8 )-( 1.12) if m E (c) = rhe a for some a > 3 or m E (c) = 0, which corresponds 
to the case "a = oo" . More precisely, we will determine radially symmetric solutions 
which converge as e — > to a solution, which does not satisfy (1.11). Moreover, for 
these solutions the discrepancy measure £ e (c £ ) does not vanish in the limit s — > 0. 

For simplicity of the following presentation we assume that z/(c) = u, p(c) = p. 
We will construct radially symmetric solutions of the form 

v(x, t) = u(r, t)e r , p(x,t) =p e (r,t), c(x,t) = c £ (r,t), fj,(x,t) = p, £ (r,t), (4.1) 



where r 



o-. If (v,p, c, p) are of this form, (1.1 )-( Jl~4 ) reduce to 



pd t u + pud T u - v-^d r (r n l d r u) 



+d r pe 

d t c £ + ud r c £ 

h 



e'^ldrC^ 2 - ed r \d r c, 



~ I 2 

= I 







m e a -^ T d r (r n 1 d r p £ ) 
-e^d r {r n - l d r c £ ) +e- 1 f'{c £ ). 



(4.2) 
(4.3) 
(4.4) 
(4.5) 
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Here we have used 

-e div ( Vc <g> Vc) 



-ediv (|<9 r c e | 2 e r 
, Ha 7, 12 







e(n — l)-\d r c £ \ e r — ed r \d r c e \ e 



since 



Ve r 



: (I-e r 



div e r 



n-1 



We note that because of (4.3) u(r,t) = ar n+1 e r for some a G R, which will be 



determined by the boundary conditions in the following. Hence we can solve (4.4) 



(4.5) together with suitable boundary conditions and c| t=0 = c 0j£ independently and 



use (4.2) afterwards to determine p e . 



4.1 Nonconvergence in the Case a = oo 



First we consider the case m = (resp. "a = oo"). In this case we consider (4.2) 



(4.5) in the domain Q = {x G M. n : \x\ > 1} together with the inflow boundary 
condition 



u(l, t) = a for all t > 0, 
c e (l,t) = l for all t > 

for some a > and the initial values 

(U,d)\ t=0 = (^T,C , e ) . 



(4.6) 
(4.7) 



Here (4.3) and (4.6) already determine u uniquely as 

a 



u(r, t) 



-n—l 



for all r > l,t> 0. 



Moreover, we choose 



Mr) = e 



r — tq 



for all r > 1 



(4.9) 



for some r > 1, where 

6 G C°°( 



such that 9(s) 



I if s < -5, 
1 if s > 5 



(4.10) 



and S G (0, ro — 1) and e G (0, 1]. Hence c is a solution of the transport equation 

d t c £ (r, t) + ^d r c E (r, t) = for r > 1, t > 0, 
c e (l,t) = l fort>0, 



4.1 Nonconvergence in the Case a = oo 
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which can be calculated with the method of characteristics. The solution for the 
initial condition above is 

-no/ ,\ ~ t ,n \coe(\ / r n - ant) if r n > ant, , t **s 
cf(r,t) := c £ (r,t) = < ' v ' ~ (4.11) 

II if r n < ant. 

By the construction we have 

/ s f-1 if r>R(t), 

s - (r - () — (i if < (4 - 12 > 

where R(t) = y/r^ + nat is the radius of the level set {c £ (x,t) = 0} = dB R ^(0). 
In order to determine p e we use that (4.2) and (4.8) imply 



d r p e = -e^\d r c £ \ 2 - ed T \d r c £ \ 2 + v d r (r n 1 d r u) - pud r i 

= -^|9 r 5 £ | 2 -^ r |9 r c £ | 2 + 9 r (t^r- 2 " +2 -^^r-"). (4.13) 

Now we decompose p £ = pi t£ + p 2 ,e + P3 such that 

d r Pi, £ (r) = -e(n - 1) ^ |<9 r c £ (r) | 2 , p 2)E (r) = -e\d r c £ (r)\ 2 for all r > 1. 
Hence up to a constant 

= a(n- l) r „ 2n+2 _ ua(n-l) 
P3 2n + 2 n 

Because of the explicit form of p £ and 

d r c £ (r,t) = ( </r n -ant-r \ ^ _ ^l-i^ (4 M) 



e 

it is easy to observe that 

p £ (r) -* £ ^ Po(r) for all r 7^ R(t) 
for some smooth p ; (1? M) \ {R(t)} — > K. Now we consider 

rR(t)+S 

[Pj,e]R(t),8 :=Pj(R{t) + 5) -pj(R(t) -6)= d r pj{s,t)ds, 

JR{t)-S 

which converges as e — » and 5 — > (in that order) to several contributions of [po] 
at For j = 2 we have that 

Mw),s = -e\dMR(t) + 5)\ 2 + e\d r c £ (R(t) - 5)\ 2 = 

if e < 5. Hence 

lim[p 2 , £ } R (t),s = 0. 



20 



4 NONCONVERGENCE 



Moreover, since p 3 is independent of e and continuous, we have 

5->0£->-0 w 



Finally, using (4.14) we obtain 

[Po]fl(t),(S 



lim[pjii(t),5 = Km[pijji(t) ) « 



n-l 

e 



pR(t)+S 


,( 


>R{t)-8 





, ( \Jr n — ant — r 



e 

2 



r ln -'\r n - ant) 



■—2 



dr 



a(n - l)r 2n - 3 (r n - ant)n- 2 
1 



r=R(t) 



a 



R(t)\ 2n ~ 2 n 



r 



R(t) 



(4.15) 



where a := J R \6'(s)\ 2 ds. Here R(t) > r for all t > and -t t -+oo oo. The 
exact solution of the classical sharp interface model, i.e., ( 1.8 )-( 1.17) with m = 
and Q£ = B ro (0) \ Si(0), v = ar 1 ""^, is given by 



v(x, t) = ar 



l-n 



X 



X 



, Q + (t) = B R{t) (0)\B 1 (0) 



where R(t) = ^jr§ + ant as before, p: Q x (0, T) — >■ R is constant in f2 ± (t) such that 

n — 1 



[p]0M) = o- 



on 0ft+(t) = Y(t). 



(4.16) 



Hence the pressure p of the limit solution as e — > differs from the solution of the 

A 1 ( R(t)\ 2n ~ 2 

sharp interface (4.16) by a time dependent factor ' — — i 



V r ° 7 



> 1, which corresponds 



to an increased surface tension coefficient that even increases strictly in time. 



Remark 4.1 From the explicit solution (|4.11|) one observes 

d r c £ (R(t),t) = 



-6'(0) 



R(t) 



r 



n—l 



Hence \d r c\ increases at the diffuse interface "r ps R(t)" as t increases, cf. Figure [T] 
Finally, we determine the limit of the discrepancy measure: 



Uc™)ipdx 



n 



M 



M 



\d r cf{r)\ 2 f{cf{r)) 



<p(r) r n dr 



f f \Jr n — ant — ro 



,n-l 



-fie 



£^0 



((TK(t) — a) 



\/ r n — ant — r$ 

e 

<p(x) dx 

R(t) 



(r n — ant) 1 n 
<p(r) r n ~ 1 dr 



(p{r) r n dr 



4.2 Nonconvergence in the Case 3 < a < oo 



21 




Figure 1: Plot of cf for t = 0, 1, 2, 3, 4 (from left to right) with a = l,e = 0.4, r 



2,n = 2. 



for all ip G C^°(f2) where a := L f(9(s)) ds, a = f R \9'(s)\ 2 ds as before, and 
0(r) = / (p(x) dx for all r e (1,M). 

JdB r (0) 

Hence 

e £ (c £ °°) ^ £ ^ - a)5 dBm in P'(^) (4.17) 

since K,(t) is strictly increasing in t > 0, we have an{t) — a ^ for all t > except 
possibly one. 



4.2 Nonconvergence in the Case 3 < a < oo 

Based on the solution for the extreme case "a = oo" from the previous section, we 
will prove essentially the same result in the case 3 < a < oo. In order to avoid 
technical difficulties with the unboundedness of {x G M. n : \x\ > 1}, we will consider 



(4.2)-(4.5) in 



n M = {x G R n : 1< \x\ < M}, 
where M > r > 1 is arbitrary, together with 

u e (l, t) = a, c e (l, t) = l for all i G (0, T), 



u e (M,t) 



M 



n— 1 



c £ (M,t) 



for all i G (0,T). 



(4.18) 
(4.19) 



Definition 4.2 (Weak Solutions) 

Let 



H (o) 



/M 
■ 



u(r)r n 1 dr — 
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be equipped with the inner product 



[u,v) H i :-- 



d r u(r)d r v(r)r n dr 



for all u, v E H^.We embed L^(1,M) > H (0) := (Hf Q) )' by identifying u E 
L 2 (1,M) with 



(0) 



[ (o)> 



M 



u(r)<f(r)r n 1 dr for all ip E H, Q y 



We call (c e ,yU e ) a weak solution of (4.4)-(4.5) together with (4. 18), (4. 19) and c £ \ t =o = 

C0,e if 

c £ - X e C([0, T]; ^(1, M)) n L 2 (0, T; F 3 (l, M)), 

cfe G L^T^J^M)), /2 e E L 2 (0,T;H\l,M)), 

where x e C°°([l, M]) unt/i = 1, x(M) = -1 and 

/•A/ pM 

{d t c £ (t),<p) H -x Hl + ar- n+1 d r c £ (r,tMr)r n - 1 dr = -m e a d^d^r"" 1 dr 

(0)' (o) J l J l 



for almost every t E (0, T) and for all ip E HLJl, M), (4.5) is satisfied pointwise 
almost everywhere, and c| i=0 = c 0>£ in H X (1,M). 

Existence of weak solutions can be proved by standard methods. E.g. it follows from 
Theorem 3.1] applied to Hi = H 1 ^, H = Hzl, 



(p(u) 



' 'e l^f + £ - 1 /o(n(r))l r^dr 



(B(«),w> 



rr-l rrl 

n (o)' n (o) 



2 

Af 



m £ a /3 / d r v(r)d r w(r)r n dr —J ad r v(r)w(r) dr 



M 



for all v,w E H\ and u E dom(<£>) := {v E H X (\,M') : ttzi L v(r)r n 1 dr = m}, 

where m := L Co j£ (r)r n ~ 1 dr, /o(s) := /(s) — f s 2 for all s G R and (3 := inf seR /"(s). 
Then / and ip are convex and the subgradient A = dip taken with respect to HZ, 
satisfies 



rr-l prl 
-°(0)'"(0) 



a+l 



Ai 



9 r (r" n+1 9 r (r"- 1 a r M))a r wr n - 1 rfr + mo£ a " 1 / /o(«(r))r n_1 dr 



A/ 



for all it G V(dip) = {v E H 3 (1,M) : j- J^ 1 v(r)r n 1 dr = m}. Then it is easy to 
verify that all conditions of [7, Theorem 3.1] are satisfied. 
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J^ 1 n £ (r)r n 1 dr in the weak formulation of the convective Cahn-Hilliard equation and 
obtain the energy identity 



Finally, if (c e , fi £ ) is a weak solution as above, we can choose (p = [i £ —[i £ with fi £ 



M 



(e\d r c £ (r,t)\ 2 + e- 1 f(c £ (r,t))) r^ 1 dr (4.20) 

+ / / m e a \V^(r,r)\ 2 r n - l drdT= / (e\d r c , £ (r)\ 2 + e" 1 f(c , £ (r))) r^ 1 dr 
Jo Ji Ji 

for all t G (0,T). 

THEOREM 4.3 L et k > 3, r G (1, M), < 5 < min(r - 1, M - r ) and T > 
such that R(T)=j/rJTnaT < M - 5, ft = {x G IR n : 1 < \x\ < M}, c , £ and 9 
be as in (4.9)-(4.10), and let (v £ ,p £ ,c £ , fi £ ) be the radially symmetric solutions of the 
form (4.1) of Ql.l[ )-(1.5), (1.7) one? boundary conditions (4.6)-(4.7) ; n- V// £ |gn = 0. 
Taen 



v e = ar 
and 

c £ -> 6 _>o 2 Xb rw (o) ~ 1 for every x G ft \ dB R{t) (0),t G (0,T), 
P.^oP m£>'(ft x (0,T)), 

where R(t) = ^/r^ + nai ; p G P'(ft x (0,T)) coincides with a function that is con- 
tinuous in x G ft \ dBn(t)(Q) for every t G (0, T), and 

[p] = K(t)o-fT on T(t) = dB m (0) for all t G (0,T) 

2n-2 



and 1 < 



V r ° 7 
^|Vc £ | 2 /(c £ 



— >t-*oo 00 ■ Moreover, 

- -^o (<r«(t) - m £>'(ft x (0, T)), (4.21) 



2 £ 

w /iere a = / R \6'(s)\ 2 ds,a = f R f(6(s))ds. 

Proof: First of all, we show that ||c £ ||j / oo((i ; m)x(o, t)) is uniformly bounded. To this 
end let W(c) be as in Section [3] Then as in (3.26) 



M 



\d r W{c e {r,t))\r n - 1 dr < C 



w ' \d r c , £ (r)\ 2 , /(co, £ (r)) , n _ x 



r"" 1 dr < C 



by (4.20) and the choice of the initial data. Hence 



sup \\W(c £ (t))\\ L ^(i,M) < sup \\d r W(c £ (t))\\ L i (hM) < C 

0<t<T,0<e<l 0<t<T,0<e<l 



24 



4 NONCONVERGENCE 



due to c £ (t, M) = 1, which implies 



SU P || c e||L°°((l,M)x(0,T)) < 
0<e<l 



(4.22) 



for some M > due to (3.27). Now let d £ := c £ — c^°, where cf is as in (4.11). First 
we will show 

>/ o 

where L 2 = L 2 ((l, M); r n_1 cfr). To this end we use that 



(d t d £ (t),<p)+m s a+l j Ad £ (r, t)A<p(r) r"" 1 dr = / &(r, t)A^(r) r n_1 dr (4.23) 
for all ip G HLs and almost every £ G (0, T), where 



and 

Moreover, 



Au(r) = r~ n+ %(r n - 1 d r u{r)) for all u G tf 2 (l, M) 
&(r,t) = m e a+1 Ac?(r,t) +m e a - 1 f'(c £ (r,t)). 



\\g e (t)\\v < m ^ +1 ||Ag-(t)|| L2 +m E a - 1 \\f (c e (t))\\ L 2 < C(T)m e a -z 
where L 2 = L 2 (l, M; r n ~ l dr) and we have used that 

\\c?(t)\\ HHl , M )<C(T)e-l 



and 



E- l \\f\c £ (t)\\h<CE £ (c £ (t))<C> 



due to \f'(s)\ 2 < C(M)f(s) for all s G [-M, M) and (4.22). Hence, choosing ip 



d £ (t) in (4.23) and integrating in time, we conclude 



sup \\d £ (t)\\ L 2 + m e 

0<t<T Jo 
rT 



" + ' I* \\Ad £ (t)f L2 dt 



pi 

< / e-^\\g £ {t)\\ L 2e^\\Ad t 
Jo 



dt. 



(4.24) 



Using the Cauchy-Schwarz and Young's inequality, we obtain 



sup \\d E (t)\\h+m e a+1 I \\Ad £ (t)\\ 2 L2 dt<C I e^" 1 !!^*)^ dt < C(T)e a - 2 



0<t<T Jo 

Combining this estimate with 

\\d r v\\ 2 L 2 < C\\v\\ L 2 \\v\\ H 2 < C\\v\\&\\Av\\v for all v G H 2 (l, M) n H%(1, M), 
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we conclude 



since a > 3. 



e\\d r d £ (t)\\l,dt < Ce^e- 1 = Ce^ 



(4.25) 



In order to determine p £ we use again (4.13) and decompose p £ = p\ t£ + p2, £ + P3 
similarly as before, where 



Ha Z |2 



<9rPi, e = -e(n - l)-|9 r q, 
Hence up to a constant 

= q(n-l) 2n+2 
F3 2n + 2 

as before. Moreover, let 

d r p? e = -e(n-l)l\d r c 



P2,e = -ed r \d r c £ \ 2 . 



vain — r 



?? 



I ; > ^oo 1 2 



PZ = -ed r \d r c 



soo|2 



be the corresponding parts of the pressure for the case a = oo. Then (4.25) implies 
that 

d r pi,e ~ d r pZ -^o in L\(1,M) x (0,T)), 

d r P2,e ~ d r p% £ in M) x (0, T)). 

Since p 3 is independent of e and the same as in the case a = oo, we conclude that 
po := limp, = limp £ °° = p™ in P'((1,M) x (0,T)), 

£—^0 e— >0 



where pf is the pressure in the case a = oo and is its limit as e — > 0. Therefore 

CPo] 



(7 



R(t)\ 2n ~ 2 n - 1 



r 



where a := L ^'(s)! 2 (is, by the result of the previous section. 



Finally, it remains to prove (4.21). First of all, because of (4.25), we conclude 
|Vc?°| 2 



lim ( e- 

e^0 



NCr 



— e- 



= in L\tt x (0,T)), 



where c^°(x,t) = c^°(|x|,t). Moreover, using (4.22) and a Taylor expansion of 
f(c £ (r,t))) around cf(r,t), we conclude 

/(c*»(r,f)) f(c e (r,t)) 



o ./l 

< 



T pM 



^1 



f{c?(r,t))d £ {r,t) 



dr dt 
drdt + C 



/•T pM 


d £ (r,t) 2 






dr dt 


/o il 







< C(M,T)e 1 ^||4||i2 ( n X ( 0)T )) + ||4||ia(nx(o,T))||/'(^ )||L 2 (nx(o ) T)) 

< C'(M,T)(e a - 3 + e^) ^ 
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since ||c £ || L oo (Cx(0T)) and \\f'{cf)\\ L 2 {nx{ ^ T)) are uniformly bounded in e G (0,1) 
due to \f (cf )\ 2 < Cf(cf). Altogether we obtain 



lim (£ E (c 

£->0 



Ucs)) = in L 1 ^ x (0,T)), 



which implies (4.21) due to (4.17). 
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